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[i]  Following  enhanced  magnetic  activity  the  fluxes  of  energetic  electrons  in  the  Earth's 
outer  radiation  belt  gradually  decay  to  quiet-time  levels.  We  use  CRRES  observations  to 
estimate  the  energetic  electron  loss  timescales  and  to  identify  the  principal  loss 
mechanisms.  Gradual  loss  of  energetic  electrons  in  the  region  3.0  <  L  <  5.0  occurs  during 
quiet  periods  ( Kp  <  3~)  following  enhanced  magnetic  activity  on  timescales  ranging  from 
1.5  to  3.5  days  for  214  keV  electrons  to  5.5  to  6.5  days  for  1.09  MeV  electrons.  The 
intervals  of  decay  are  associated  with  large  average  values  of  the  ratio  fpJfce  (>7), 
indicating  that  the  decay  takes  place  in  the  plasmasphere.  We  compute  loss  timescales  for 
pitch-angle  scattering  by  plasmaspheric  hiss  using  the  PADIE  code  with  wave  properties 
based  on  CRRES  observations.  The  resulting  timescales  suggest  that  pitch  angle  scattering 
by  plasmaspheric  hiss  propagating  at  small  or  intermediate  wave  normal  angles  is 
responsible  for  electron  loss  over  a  wide  range  of  energies  and  L  shells.  The  region 
where  hiss  dominates  loss  is  energy-dependent,  ranging  from  3.5  <  L  <  5.0  at  214  keV  to 
3.0  <  L  <  4.0  at  1.09  MeV.  Plasmaspheric  hiss  at  large  wave  normal  angles  does  not 
contribute  significantly  to  the  loss  rates.  At  E  =  1.09  MeV  the  loss  timescales  are 
overestimated  by  a  factor  of  ~5  for  4.5  <  L  <  5.0.  We  suggest  that  resonant  wave-particle 
interactions  with  EMIC  waves,  which  become  important  at  MeV  energies  for  larger 
L  (L  >  ~4.5),  may  play  a  significant  role  in  this  region. 
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1.  Introduction 

[2]  Energetic  electrons  (E  >  100  keV)  in  the  Earth’s 
radiation  belts  are  generally  confined  to  two  distinct 
regions.  The  inner  radiation  belt  lies  in  the  range  1.2  <  L 
<  2.5  and  exhibits  long-term  stability.  In  contrast,  the  outer 
radiation  belt,  which  lies  in  the  range  3  <  L  <  7,  is  highly 
dynamic,  particularly  during  enhanced  magnetic  activity 
[e.g.,  Paulikas  and  Blake ,  1979;  Baker  et  al. ,  1986,  1994, 
1997;  Li  et  al .,  1997;  Reeves  et  al .,  1998].  This  variability  is 
caused  by  an  imbalance  between  acceleration  and  loss 
processes  both  of  which  tend  to  be  enhanced  during 
magnetically  disturbed  periods  [e.g.,  Summers  et  al ., 
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2004].  Understanding  this  variability  is  important  since 
enhanced  fluxes  of  relativistic  electrons  (£  >  1  MeV) 
damage  satellites  in  Earth  orbit  and  pose  a  risk  to  humans 
in  space.  Indeed,  enhanced  fluxes  of  relativistic  electrons 
have  been  associated  with  a  number  of  spacecraft  anomalies 
and  even  failures  [Baker  et  al .,  1998a;  Baker ,  2001]. 
Consequently,  predicting  their  appearance  has  become  one 
of  the  outstanding  challenges  of  magnetospheric  physics. 
Furthermore,  energetic  electrons  can  penetrate  to  low  alti¬ 
tudes  where  they  affect  the  ionization,  conductivity,  and 
chemistry  of  the  stratosphere  and  mesosphere  [e.g.,  Thorne , 
1980;  Lastovicka ,  1996;  Callis  et  al .,  1998],  thereby  pro¬ 
viding  an  important  coupling  mechanism  between  the 
magnetosphere  and  the  middle  atmosphere. 

[3]  Electrons  with  energies  up  to  a  few  hundred  keV  are 
injected  into  the  inner  magnetosphere  during  storms  and 
substorms.  Injection  of  electrons  in  the  energy  range  10 
100  keV  into  the  outer  zone  leads  to  the  excitation  of  intense 
whistler  mode  chorus  waves  outside  the  plasmasphere  on 
the  dawnside  of  the  magnetosphere  [e.g.,  Tsurutani  and 
Smith ,  1977;  Meredith  et  al .,  2001,  2003a].  At  higher 
energies,  injected  electrons  with  energies  in  the  range 
100-300  keV  form  a  seed  population  of  electrons  [e.g., 
Baker  et  al .,  1998b;  Ohara  et  al .,  2000]  which  may 
subsequently  be  accelerated  to  relativistic  energies  by  pro¬ 
cesses  acting  within  the  magnetosphere  itself  [Li  et  al ., 
1997].  While  several  acceleration  mechanisms  have  been 
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proposed  (see  the  reviews  by  Li  and  Temerin  [2001],  Horne 
[2002],  and  Friedel  et  al.  [2002]),  recent  experimental  and 
modelling  work  suggests  that  the  acceleration  is  caused  by  a 
combination  of  enhanced  inward  radial  diffusion  driven  by 
ULF  waves  [e.g.,  Elkington  et  al .,  1999;  Liu  et  al .,  1999; 
Hudson  et  al.,  2000]  and  local,  chorus-driven  acceleration 
outside  the  plasmapause  [e.g.,  O  'Brien  et  al,  2003;  Horne 
et  al.,  2006;  Summers  et  al .,  2005;  Varotsou  et  al.,  2005]. 
However,  during  quiet  conditions  following  enhanced  mag¬ 
netic  activity,  the  acceleration  processes  are  ineffective  and 
the  energetic  electrons  gradually  decay  to  their  quiet  time 
values.  In  order  to  understand  radiation  belt  variability  and 
develop  more  realistic  physical  models,  it  is  essential  to 
understand  and  quantify  the  decay  processes.  The  periods  of 
gradual  decay  provide  a  unique  possibility  to  isolate  the  loss 
processes  acting  during  these  intervals  and  to  determine  the 
decay  lifetimes. 

[4]  There  are  several  loss  mechanisms  that  operate  during 
quiet  time  decay.  Coulomb  collisions  with  atmospheric 
constituents  are  important  in  the  inner  zone  and  are  the 
dominant  loss  mechanism  closest  to  the  Earth,  while  reso¬ 
nant  interactions  with  plasma  waves  become  increasingly 
important  farther  out.  The  region  where  Coulomb  collisions 
dominate  is  energy-dependent,  ranging  from  L  <  2.2  for 
100  keV  electrons  to  L  <  1.4  for  1.5  MeV  electrons  [Abel 
and  Thorne ,  1998].  Ground-based  VLF  transmitters,  used 
for  communications  with  submarines,  leak  out  of  the 
atmosphere  at  night  and  propagate  in  the  inner  magneto¬ 
sphere  in  the  whistler  mode.  These  waves  are  most 
important  just  outside  the  region  dominated  by  Coulomb 
collisions.  Using  average  wave  characteristics  as  input  to 
their  model,  Abel  and  Thorne  [1998]  showed  that  pitch- 
angle  scattering  by  VLF  transmitters  is  likely  to  be  the  most 
effective  loss  mechanism  for  100  keV  electrons  in  the 
region  from  2.2  <  L  <  2.8  and  for  1.5  MeV  electrons  in 
the  region  1 .4  <  L  <  2.2.  Lightning-generated  whistlers  are 
most  effective  just  outside  the  region  dominated  by  VLF 
transmitters,  being  most  effective  from  2.8  <  L  <  4.4  for  100 
keV  electrons  and  from  2.2  <  L  <  2.6  for  1.5  MeV  electrons 
[Abel  and  Thorne ,  1998].  Farther  out,  but  remaining  inside 
the  plasmapause,  resonant  interactions  with  both  plasma- 
spheric  hiss  and  electromagnetic  ion  cyclotron  (EMIC) 
waves  are  believed  to  dominate  [Thorne  et  al.,  1973;  Lyons 
et  al.,  1972;  Albert,  1994;  Summers  and  Thorne,  2003; 
Summers,  2005].  Outside  of  the  plasmapause,  whistler 
mode  chorus  waves  can  contribute  effectively  to  the  loss 
of  energetic  electrons  [Horne  and  Thorne,  2003;  O'Brien  et 
al.,  2004;  Thorne  et  al.,  2005a;  Summers  et  al.,  2005]. 

[5]  The  studies  described  above  have  identified  loss 
timescales  in  various  regions  but  have  been  limited  in  the 
coverage  of  energies  and  L  shells.  Here  we  provide  a  more 
extensive  survey  using  the  CRRES  data  set.  The  objectives 
are  to  determine  experimentally  the  loss  timescales  for 
energetic  electrons  over  a  wide  range  of  energies  and  L 
shells,  to  determine  the  conditions  associated  with  loss,  and 
to  identify  the  principal  loss  mechanisms. 

2.  Instrumentation 

[6]  The  Combined  Release  and  Radiation  Effects  Satellite 
(CRRES)  is  particularly  well-suited  to  studies  of  wave- 
particle  interactions  in  the  radiation  belts  both  because  of 


its  orbit  and  sophisticated  suite  of  wave  and  particle  instru¬ 
ments  [Johnson  and  Kierein,  1992].  The  spacecraft  was 
launched  on  25  July  1990  and  operated  in  a  highly  elliptical 
geosynchronous  transfer  orbit  with  a  perigee  of  305  km,  an 
apogee  of  35,768  km,  and  an  inclination  of  18°.  The  orbital 
period  was  approximately  10  hours,  and  the  initial  apogee 
was  at  a  magnetic  local  time  (MLT)  of  0800  MLT.  The 
magnetic  local  time  of  apogee  decreased  at  a  rate  of 
approximately  1.3  hours  per  month  until  the  satellite  failed 
on  1 1  October  1991,  when  its  apogee  was  at  about  1 400  MLT. 
The  satellite  swept  through  the  heart  of  the  radiation  belts  on 
average  approximately  5  times  per  day,  providing  good 
coverage  of  this  important  region  for  almost  15  months. 

[7]  The  electron  data  used  in  this  study  were  collected  by 
the  Medium  Electrons  A  (MEA)  experiment.  This  instru¬ 
ment,  which  used  momentum  analysis  in  a  solenoidal  field, 
had  17  energy  channels  ranging  from  153  keV  to  1.58  MeV 
[Vampola  et  al.,  1992].  The  wave  data  used  in  this  study 
were  provided  by  the  Plasma  Wave  Experiment.  This 
experiment  provided  measurements  of  electric  fields  from 
5.6  Hz  to  400  kHz,  using  a  100  m  tip-to-tip  long  wire 
antenna,  with  a  dynamic  range  covering  a  factor  of  at  least 
105  in  amplitude  [Anderson  et  al.,  1992].  The  electric  field 
detector  was  thus  able  to  detect  waves  from  below  the  lower 
hybrid  resonance  frequency  {[lur)  to  well  above  the  upper 
hybrid  resonance  frequency  {fuHR)  for  a  large  fraction  of 
each  orbit. 

3.  Data  Analysis 

3.1.  CRRES  Database 

[8]  In  order  to  study  electron  loss  rates  and  the  role  of 
plasmaspheric  hiss  in  the  loss  process,  we  constructed  a 
database  from  the  wave  and  particle  data.  The  electron 
differential  number  flux  perpendicular  to  the  ambient  mag¬ 
netic  field  for  each  energy  level  of  the  MEA  instrument  and 
the  ratio  fpJfce,  together  with  the  magnetic  field  intensities 
in  the  range  0.1  </<  2  kHz  and  the  electric  field  intensities 
in  the  range  fce  <  f<  2 f  .e  were  rebinned  as  a  function  of  half 
orbit  (outbound  and  inbound)  and  L  in  steps  of  0. 1  L  as 
detailed  by  Meredith  et  al.  [2004].  The  data  were  recorded 
together  with  the  universal  time  (UT),  magnetic  latitude 
(Xm),  magnetic  local  time  (MLT),  and  time  spent  in  each  bin 
with  the  same  resolution.  The  subsequent  analysis  of  the 
particle  data  was  restricted  to  the  equatorial  perpendicular 
fluxes,  defined  to  be  observations  within  ±15°  of  the 
magnetic  equator  to  reduce  latitudinal  effects  in  measure¬ 
ments  of  the  perpendicular  flux  and  the  ratio  fpJf.e.  If  we 
assume  a  dipole  field  then  this  criterion  restricts  the  analysis 
to  equatorial  pitch  angles,  a,  in  the  range  (60°  <  aeo  < 
120°). 

3.2.  Determination  of  the  Loss  Timescale 

[9]  During  quiet  periods  the  energetic  electron  flux,  J, 
often  exhibits  an  exponential  decay.  To  quantify  the  decay 
and  compare  results  for  different  energies  and  L  shells,  we 
have  calculated  a  loss  timescale,  t,  by  fitting  an  exponential 
function  of  the  form  J  =  Ae~*//T)  to  periods  of  gradual  decay. 
To  prevent  the  fits  from  being  dominated  by  the  most 
intense  fluxes,  we  obtain  the  decay  time  constant  by  fitting 
a  linear  function  to  the  natural  logarithm  of  the  flux.  The 
fitting  intervals  are  selected  automatically  to  prevent  ob- 
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server  bias  from  influencing  the  results.  We  measure  the 
strength  of  the  linear  relationship  between  the  natural 
logarithm  of  the  flux  and  the  decay  time  using  the  Pearson 
correlation  coefficient.  The  value  of  this  coefficient  varies 
from  -1  to  +1  with  —1  indicating  a  perfect  negative 
(inverse)  correlation,  0  no  correlation,  and  +1  a  perfect 
positive  correlation.  We  proceed  as  follows.  The  Pearson 
correlation  coefficient  is  determined  for  the  first  ten  points. 
If  the  Pearson  coefficient  is  negative  with  an  absolute  value 
less  than  0.95  or  positive  the  fit  is  discarded  and  the  starting 
point  is  incremented  by  one  point.  The  process  is  then 
repeated  throughout  the  data  set.  When  the  Pearson  corre¬ 
lation  coefficient  is  negative  with  an  absolute  value  greater 
than  0.95  the  number  of  points  included  in  the  fit  is 
increased  in  unit  increments  up  to  a  maximum  of  75  points. 
The  fitting  interval  is  then  chosen  to  be  the  fit  with  the 
largest  absolute  value  of  the  Pearson  coefficient.  The  whole 
process  is  then  repeated  starting  from  the  next  point  after  the 
fitting  interval.  The  process  is  repeated  for  the  whole  data 
set  as  a  function  of  time  for  selected  energies  and  L  shells. 

[io]  We  tested  the  sensitivity  of  the  technique  to  the 
minimum  number  of  points  in  the  fit  by  using  a  minimum 
of  five  points  as  compared  with  ten  and  found  that  the 
results  were  relatively  insensitive  to  the  minimum  number 
of  point  in  the  fit.  The  nature  of  the  fitting  technique  is  such 
that  we  do  not  fit  to  intervals  with  less  than  five  points 
which  corresponds  to  a  minimum  duration  of  approximately 
25  hours  due  to  the  orbital  period  of  CRRES.  Periods  of  fast 
decay,  associated  with  the  main  phase  of  geomagnetic 
storms,  are  thus  excluded  from  the  fits  and  we  concentrate 
on  the  periods  of  more  gradual  decay  following  enhanced 
magnetic  activity. 

3.3.  Determination  of  the  Plasma  Environment 

[n]  The  ratio  of  the  electron  plasma  frequency,^,  to  the 
electron  gyro  frequency,  f.e ,  is  an  important  parameter  for 
electron  acceleration  and  loss  in  the  inner  magnetosphere 
[Summers  et  al.,  1998;  Horne ,  2002].  Given  the  presence  of 
the  appropriate  whistler  mode  waves,  low  values  of  this 
ratio  are  associated  with  the  plasma  trough  and  energy  and 
pitch-angle  diffusion  [Summers  et  al .,  1998;  Horne  et  al ., 
2003].  High  values  of  this  ratio  are  associated  with  the 
plasmasphere,  pitch-angle  scattering  and  subsequent  loss  to 
the  atmosphere  [Lyons  et  al .,  1972;  Summers  et  al .,  1998; 
Summers ,  2005]. 

[12]  Loss  processes  inside  the  plasmapause  may  be  very 
different  to  those  outside  and  have  different  relative  impor¬ 
tance.  It  is  therefore  helpful  to  identify  the  location  of  the 
observations  with  respect  to  the  plasmapause,  which  can  be 
inferred  from  measurements  of  the  ratio  o  ffpJfce.  In  order  to 
construct  their  empirical  plasmaspheric  density  model, 
Sheeley  et  al.  [2001]  adopted  an  L-shell  dependent  bound¬ 
ary  number  density,  given  by  nh  =  107(6.6/£)4  m“3,  to 
separate  plasmaspheric-like  and  trough-like  values.  Densi¬ 
ties  below  nh  were  considered  trough-like  and  densities 
above  nh  were  considered  to  be  plasmaspheric-like.  We  use 
the  same  criterion  and  express  the  boundary  ratio  ifpJfce)h  in 
terms  of  the  boundary  number  density,  nh,  and  the  local 
magnetic  field  strength,  B,  as 

(fpe/fce)b=  \/(nhme/e.0B2)  (1) 


If  we  assume  a  dipole  field,  then 

B  =  Beq  yj(  I  +  3cos20)/L3  sin6  0  (2) 

where  Beq  is  the  mean  value  of  the  magnetic  field  on  the 
equator  at  the  surface  of  the  Earth,  and  0  is  the  colatitude. 
Substituting  in  values  for  the  appropriate  constants,  and 
keeping  the  plasma  density  constant  with  magnetic  latitude, 
we  obtain 

(ft*/fce)h=  1-416Z,  sin6  e/\/(l  +  3  cos2  0)  (3) 

The  boundary  value,  separating  the  plasmasphere-like 
values  from  trough-like  values,  for  measurements  made 
within  15°  of  the  magnetic  equator  on  a  given  geomagnetic 
field  line  thus  lies  in  the  range  1.23L  ±  0.1 8L.  Typically,  in 
the  plasmasphere,  the  ratio  fpJfce  is  much  larger  than  the 
boundary  value,  and,  in  the  plasma  trough,  the  ratio  fpJfce  is 
much  smaller  than  the  boundary  value. 

4.  Results 
4.1.  £  =  3.55 

[13]  The  Kp  index  and  the  equatorial  perpendicular  dif¬ 
ferential  number  flux  of  1.09  MeV  electrons  at  L  =  3.55  are 
plotted  as  a  function  of  day  number  from  01/01/90  for  the 
entire  CRRES  mission  in  Figure  1.  The  Kp  index  has  been 
smoothed  using  a  15  hour  running  mean  and  the  electron 
data  have  been  color-coded  according  to  the  value  of  the 
ratio  fpJfce-  In  the  equatorial  region  at  L  =  3.55,  (fptJfce)h  = 
4.4  ±  0.6  so  that  measurements  made  when  fpJfce  is  above 
or  below  ~4.4  may  be  regarded  as  plasmaphere-like  or 
trough-like,  respectively. 

[14]  Large  increases  in  the  flux  of  relativistic  electrons  are 
associated  with  increased  magnetic  activity  as  monitored  by 
the  Kp  index,  and,  in  particular,  when  Kp  >  5.  These  events 
are  associated  with  magnetic  storms,  as  monitored  by  the 
Dst  index.  The  increase  can  be  very  rapid,  such  as  the  event 
on  day  number  332,  or  may  take  place  over  a  timescale  of 
days.  The  very  rapid  events  may  be  attributed  to  either 
shock  acceleration  [Li  et  al.,  1993;  Hudson  et  al.,  1997]  or 
the  inward  movement  of  the  inner  edge  of  the  outer  belt. 
The  longer-lasting  flux  increases  take  place  over  a  period  of 
2-4  days  and  are  associated  with  low  values  offpe/f  e  and 
hence  the  trough  region.  These  events  are  known  to  be 
associated  with  enhanced  levels  of  magnetic  activity,  en¬ 
hanced  chorus  activity,  and  enhanced  fluxes  of  seed  elec¬ 
trons  [e.g.,  Meredith  et  al.,  2002a,  2002b,  2003b;  Miyoshi  et 
al.,  2003]. 

[15]  In  the  absence  of  further  significant  magnetic  activity 
(Kp  <  ~3)  the  elevated  poststorm  relativistic  electron  fluxes 
gradually  decay  to  quiet  time  values.  During  quiet  periods 
the  fluxes  may  fall  by  as  much  as  two  orders  of  magnitude 
over  a  period  of  ~20  days.  The  best  fits  to  the  selected 
decay  periods  are  overplotted  on  the  data  and  yield  a  decay 
time  constant  of  5.7  ±  0.6  days,  where  the  quoted  error,  here 
and  henceforth,  is  the  error  of  the  mean.  These  periods  of 
gradual  decay  are  all  associated  with  elevated  values  of  the 
ratio  fpJfce.  The  average  value  of  fpJfce  for  the  fitting 
intervals  is  9.5  ±  0.2,  indicating  that  these  periods  of  gradual 
flux  decay  take  place  largely  in  the  plasmasphere.  The 
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Kp  index  and  the  equatorial  1.090  MeV  flux  at  L  =  3.55 
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Figure  1.  The  Kp  index  and  the  equatorial  perpendicular  differential  number  flux  of  1 .09  MeV  electrons 
at  L  =  3.55  as  a  function  of  day  number  for  the  entire  CRRES  mission.  The  Kp  index  has  been  smoothed 
using  a  15-hour  running  mean  and  the  data  points  are  color-coded  according  to  the  value  o ffpe/fce  at  the 
time  of  the  measurement.  The  Kp  value  of  3“  is  represented  by  the  horizontal  dashed  line.  The  best  fits  to 
the  selected  decay  periods  are  overplotted  on  the  data. 


average  value  of  the  Kp  index  during  the  intervals  of  decay 
is  2.23  ±  0.02,  indicative  of  relatively  quiet  magnetic 
activity. 

4.2.  L  -  4.55 

[16]  The  Kp  index  and  the  equatorial  perpendicular  differ¬ 
ential  number  flux  of  1.09  MeV  electrons  at  L  =  4.55 
are  plotted  as  a  function  of  day  number  from  01/01/90  in 
Figure  2.  In  the  equatorial  region  at  L  =  4.55  measurements 
made  when  fpJfce  is  above  or  below  ~5.6  may  be  regarded  as 
plasmaphere-like  or  trough-like,  respectively.  At  L-  4.55  we 
see  a  larger  number  of  flux  enhancements  due  to  the  fact  that 
weaker  storms  can  result  in  flux  increases  at  larger  L.  The 
increases  are  again  associated  with  low  values  of  fpJfce , 
typically  fpJfce  <  4.  At  L  =  4.55  the  decay  time-constant  is 
5.3  ±  0.6  days,  comparable  to  the  timescale  at  L  =  3.55.  The 
average  value  o ifpJfce  during  the  decay  periods  is  1 1 .0  ±  0.4, 
consistent  with  electron  decay  in  this  region  taking  place  in 
the  plasmasphere.  The  average  value  of  the  Kp  index  during 
the  intervals  of  decay  is  1 .83  ±  0.02,  indicative  of  relatively 
quiet  magnetic  activity. 

[17]  The  Kp  index  and  the  equatorial  perpendicular  elec¬ 
tron  differential  number  flux  of  214  keV  electrons  at  L  = 
4.55  are  shown  for  comparison  with  the  MeV  electrons  in 
Figure  3.  During  high  magnetic  activity  the  ratio  fpJfce  is 
low  denoting  the  plasma  trough.  The  flux  levels  also 
tend  remain  high.  The  loss  timescales  are  more  rapid  than 
at  1.09  MeV,  and,  during  quiet  periods,  the  fluxes  may  fall 
by  as  much  as  two  orders  of  magnitude  over  a  period  of 
~6  days.  The  best  fits  to  the  selected  decay  periods  yield  a 
decay  time  constant  of  2.0  ±  0. 1  days,  a  factor  of  2.65  times 


faster  than  at  1 .09  MeV.  The  average  value  of  the  ratio  fpJfce 
for  the  fitting  intervals  is  9.9  ±  0.3,  consistent  with  electron 
decay  at  this  energy  also  occurring  in  the  plasmasphere.  The 
intervals  of  decay  again  occur  during  quiet  magnetic  con¬ 
ditions  when  the  average  value  of  the  Kp  index  is  1 .7  ±  0. 1 . 

4.3.  3  <  L  <  5 

[is]  The  analysis  was  repeated  for  different  energies  and 
L  shells  and  the  resulting  decay  timescales  are  plotted  as  a 
function  of  energy  and  L  shell  for  the  region  3  <  L  <  5  in 
Figure  4.  At  1.09  MeV  the  decay  timescales  lie  in  the  range 
5.5 -6.5  days  and  show  little  variation  with  L  shell.  In 
contrast,  at  510  keV  the  decay  timescales  are  less,  ranging 
from  2.5  to  5.5  days,  and  display  a  general  tendency  to 
increase  with  increasing  L  for  L  >  3.6.  The  shortest  time- 
scales  are  seen  at  2 14  keV,  ranging  from  1 .5  to  3.5  days,  and 
display  a  tendency  to  increase  with  increasing  L  for  L  >  4. 
Beyond  L  =  ~5  the  data  exhibit  more  variability  and  it  is  not 
possible  to  calculate  reliable  average  loss  timescales  using 
our  fitting  criteria.  The  fitting  technique  also  breaks  down 
inside  L  =  3.3  for  2 14  keV  electrons,  due  to  increased  scatter 
between  neighboring  data  points. 

[i9]  The  average  value  of  the  ratio  fpJfce  during  the 
selected  intervals  of  gradual  decay  for  214  keV,  510  keV, 
and  1 .09  MeV  electrons  at  each  L  shell  studied  is  shown  in 
Figure  5.  The  dashed  line  represents  the  boundary  value 
between  the  trough  and  plasmaspheric-like  material  with 
trough-like  material  below  the  line  and  plasmaspheric-like 
material  above  the  line.  Here  fpJfce  lies  well  above  the  line 
for  all  energies  and  L  shells  suggesting  that  the  intervals  of 
gradual  decay  take  place  in  the  plasmasphere. 
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Figure  2.  The  Kp  index  and  the  equatorial  perpendicular  differential  number  flux  of  1 .09  MeV  electrons 
at  L  =  4.55  in  the  same  format  as  Figure  1. 


[20 J  The  average  value  of  the  Kp  index  during  the 
selected  intervals  of  gradual  decay  for  214  keV,  510  keV, 
and  1 .09  MeV  electrons  at  each  L  shell  studied  is  shown  in 
Figure  6.  The  intervals  of  decay  are  associated  with  low 
values  of  Kp  (Kp  <  3~),  indicative  of  quiet  time  conditions. 
There  is  also  a  tendency  for  the  average  value  of  Kp  during 


the  decay  periods  to  decrease  with  increasing  L  shell.  Since 
the  quiet-time  decay  takes  place  largely  in  the  plasma- 
sphere,  quiet-time  decay  only  becomes  possible  at  a  given 
location  once  the  plasmapause  has  expanded  beyond  the 
location,  which  can  only  occur  when  Kp  becomes  small 
[' Carpenter  and  Anderson ,  1992].  On  average,  quieter 
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Figure  3.  The  Kp  index  and  the  equatorial  perpendicular  differential  number  flux  of  214  keV  electrons 
at  L  =  4.55  in  the  same  format  as  Figure  1. 
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Figure  4. 

electrons. 


Loss  timescole  versus  L  shell 


Measured  electron  loss  timescales  versus  L  shell  for  214  keV,  510  keV,  and  1.09  MeV 


Ratio  fpe/fce  versus  L  shell 


Figure  5.  The  average  value  of  fpJfce  during  the  selected  intervals  of  decay  versus  L  shell  for  2 14  keV, 
5 1 0  keV,  and  1 .09  MeV  electrons.  The  dashed  line  represents  the  boundary  value  between  trough  and 
plasmaspheric-like  material. 
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Figure  6.  The  average  value  of  the  Kp  index  during  the  selected  intervals  of  decay  versus  L  shell  for 
214  keV,  510  keV,  and  1.09  MeV  electrons. 

magnetic  conditions  and  smaller  values  of  Kp  are  required 
with  increasing  L  which  explains  the  observed  trend. 

5.  Loss  Mechanisms 

[21]  The  observations  of  gradual  periods  of  electron  loss 
lasting  for  several  days  or  more  are  largely  associated  with 
high  values  of fpJfce  and  hence  take  place  in  the  plasma- 
sphere.  While  observations  on  the  duskside  cannot  confirm 
that  the  electrons  spend  all  of  their  drift  orbits  inside  the 
plasmapause,  the  observations  of  loss  associated  with  large 
values  of  fpJf  e  on  the  dawnside  suggest  that  during  the 
periods  of  gradual  decay,  the  electrons  do  indeed  remain 
inside  the  plasmapause  for  the  majority  of  their  drift  orbits. 

[22]  The  chorus  waves,  thought  to  play  a  role  in  the 
acceleration  of  electrons  to  relativistic  energies  in  the  outer 
radiation  belt,  tend  to  be  confined  to  the  low-density  regions 
outside  of  the  plasmapause  and  so  are  unlikely  to  play  a  role 
in  electron  loss  in  the  plasmasphere.  On  the  other  hand, 
electromagnetic  ion  cyclotron  (EMIC)  waves,  plasma- 
spheric  hiss,  lightning-generated  whistlers,  and  VLF  trans¬ 
mitter  signals  are  observed  in  the  plasmasphere  and  may  all 
contribute  to  pitch-angle  scattering  and  subsequent  electron 
loss  in  this  region  [e.g.,  Albert ,  2003;  Summers  and  Thorne , 

2003;  Thorne  et  al. ,  2005b;  Summers,  2005]. 

[23]  Electromagnetic  ion  cyclotron  waves,  which  propa¬ 
gate  in  bands  below  the  proton  gyrofrequency,  are  able  to 
interact  with  energetic  electrons  resulting  in  pitch-angle 
scattering  and  loss  to  the  atmosphere  [e.g.,  Summers  and 
Thorne,  2003;  Albert,  2003].  However,  electron  minimum 
resonant  energies  are  only  observed  to  fall  below  2  MeV  in 
the  region  L  >  4.5  when  fpJfce  >  10  and  largely  occur  on  the 


duskside  [Meredith  et  al.,  2003c].  Furthermore,  these  lower- 
energy  scattering  events  tend  to  be  associated  with  magnetic 
storms.  It  is  therefore  unlikely  that  EMIC  waves  will  play 
the  most  significant  role  in  the  loss  of  electrons  with 
energies  in  the  range  100  keV  <  E  <  1  MeV  during  the 
relatively  quiet  periods  associated  with  gradual  electron 
decay  and,  in  particular,  inside  of  L  =  4.5. 

[24]  Plasmaspheric  hiss  is  a  broadband,  structureless 
emission  which  occurs  in  the  frequency  range  from 
100  Hz  to  several  kHz.  This  whistler  mode  emission  is 
observed  inside  of  the  plasmapause  and  is  most  intense 
during  storms,  although  the  emissions  also  persist  during 
relatively  quiet  times  [e.g..  Smith  et  al.,  1974;  Thorne  et  al., 
1977;  Meredith  et  al.,  2004].  The  minimum  resonant 
energies  for  a  band  of  hiss  in  the  range  3.0  <  L  <  5.0  are 
typically  less  than  100  keV,  and  thus  the  waves  can  resonate 
with  electrons  with  E  >  100  keV.  Plasmaspheric  hiss  could 
thus  play  an  important  role  in  electron  loss  inside  the 
plasmapause  and,  in  particular,  to  the  gradual  loss  observed 
when  quiet  periods  follow  geomagnetic  storms. 

6.  Estimation  of  Loss  Rates  Due  to 
Plasmaspheric  Hiss 

[25]  We  investigate  the  role  of  plasmaspheric  hiss  as  a 
loss  process  using  wave  observations  from  the  CRRES 
spacecraft  to  calculate  pitch-angle  diffusion  rates  for  elec¬ 
trons.  The  diffusion  rates  are  calculated  using  the  PADIE 
(Pitch  Angle  and  Energy  Diffusion  of  Ions  and  Electrons) 
code  [Glauert  and  Horne,  2005]. 

[26]  Since  resonant  scattering  by  hiss  is  not  sensitive  to 
the  ion  composition,  an  electron/proton  plasma  is  assumed. 


Kp  index  versus  L  shell 


7  of  13 


A05212 


MEREDITH  ET  AL.:  ENERGETIC  ELECTRON  LOSS  TIMESCALES 


A05212 


The  determination  of  the  diffusion  coefficients  then  requires 
knowledge  of  the  distribution  of  the  wave  power  spectral 
density  with  frequency  and  wave  normal  angle,  together 
with  the  ratio  fpJfce,  wave  mode,  and  the  number  of 
resonances.  We  calculate  the  bounce-averaged  pitch-angle 
diffusion  coefficients  for  whistler  mode  hiss  for  Landau 
(n  =  0)  and  ±5  cyclotron  harmonic  resonances. 

[27]  The  waves  are  assumed  to  have  a  Gaussian  frequency 
distribution  given  by 


[  A2  exp^- 

-w: 

)  U )ic  <U<  uw 

(4) 

lo 

otherwise, 

equation  (7)  becomes 

#  f/n  v-r  0  dg\  T sin 2uo 

(in 

Equation  (11),  together 

with  the  boundary  conditions 

g(nL)  =  0. 

(12) 

da 0  V2/ 

(13) 

where  B 2  is  the  power  spectral  density  of  wave  magnetic 
field  (in  T2  Hz-1),  u;m  and  are  the  frequency  of 
maximum  wave  power  and  bandwidth,  respectively,  uj/c.  and 
ujmc  are  lower  and  upper  bounds  to  the  wave  spectrum 
outside  which  the  wave  power  is  zero,  and  A2  is  a 
normalization  constant  given  by 


a>JM 


bu)  y/TV 


2  /  g sina0</a0  =  1,  (14) 

Jc\L 

where  aL  is  the  loss  cone  pitch-angle,  can  be  cast  as  a  two- 
point  boundary  value  problem  in  four  variables  [ Albert , 
1994],  which  can  be  solved  to  obtain  the  lifetime,  t.  We 
solved  this  boundary  value  problem  with  a  routine  from  the 
NAG  library,  D02GAF,  which  uses  a  finite  difference 
method  combined  with  a  deferred  correction  technique  and 
a  Newton  iteration. 


where  Z?H  is  the  wave  amplitude  in  units  of  Tesla.  The 
distribution  of  wave  normal  angles  i|;  is  also  assumed  to  be 
Gaussian,  given  by 


f'xp( 

-(  «•*)) 

(6) 

[0 

otherwise, 

where  X  =  tan(ij>),  bX  is  the  angular  width,  X„,  is  the  peak, 
and  Xlc  and  Xuc  are  the  lower  and  upper  bounds  to  the  wave 
normal  distribution  outside  of  which  the  wave  power  is 
zero. 

[28]  Once  the  pitch-angle  diffusion  rate  is  calculated,  the 
timescale  for  the  electrons  to  pitch-angle  scatter  into  the  loss 
cone  can  be  determined.  Following  previous  work  [Lyons  et 
al.,  1972;  Albert,  1994],  we  assume  that  the  electron 
distribution  function,  f,  satisfies  the  one-dimensional 
pitch-angle  diffusion  equation, 

where  is  the  equatorial  pitch-angle,  (Dtm(oto))  is  the 
bounce-averaged  pitch-angle  diffusion  coefficient  in  units 
of  s-1,  and 


r(oo)  =  1.30  -  0.56  sin  oto  (8) 


is  an  approximation  to  the  mirror  latitude  dependence  of  the 
bounce  period. 

[29]  By  assuming  that  /  can  be  factorized  into  time- 
dependent  and  pitch-angle  dependent  functions, 

/(<W)  =  F(t)g(  oo),  (9) 

and  that  the  precipitation  lifetime,  t  is  given  by 


-F 

dF/dt ' 


(10) 


6.1.  Wave  Model 

[30]  We  use  a  model  based  on  CRRES  observations  for 
the  plasmaspheric  hiss  wave  intensities.  Since  whistler 
mode  chorus  waves  can  fall  into  the  hiss  band  outside  of 
L  ~  3.5,  we  adopt  a  criterion  based  on  the  wave  amplitude 
in  the  ECH  band  (fcc  <  f  <  2 fce)  to  distinguish  between 
plasmaspheric  hiss  and  chorus.  Specifically,  we  use  the 
criterion  that  the  ECH  wave  amplitude  for  frequencies  in 
the  range  fce  <  f  <  2 \fce  must  be  less  than  0.0005  mV  itT  1  in 
order  for  wave  emissions  below  fce  in  the  frequency  range 
0.1  <  f  <  2  kHz  to  be  identified  as  plasmaspheric  hiss 
[Meredith  et  al.,  2004]. 

[31]  Energetic  electrons  ( E  >  200  keV)  in  the  Earth’s 
outer  radiation  belt  drift  around  the  Earth  on  timescales  of 
the  order  of  1  hour  or  less  which  means  that  they  typically 
complete  many  orbits  during  the  periods  of  gradual  decay. 
The  electrons  interact  with  plasmaspheric  hiss  while  their 
drift  orbit  lies  inside  the  plasmapause  which,  during  quiet 
conditions  in  the  region  3  <  L  <  5,  is  typically  the  entire 
orbit.  Therefore  we  require  a  global  model  of  the  hiss 
intensities  to  obtain  an  estimate  of  the  loss  rates.  Hiss 
intensities  vary  as  a  function  of  spatial  location  and  mag¬ 
netic  activity  [Meredith  et  al.,  2004].  Since  the  intervals  of 
gradual  decay  take  place  during  quiet  periods  ( Kp  <  3  ),  the 
survey  is  restricted  to  measurements  taken  when  the  instan¬ 
taneous  value  of  Kp  is  less  than  3~.  The  average  value  of 
the  intensity  of  plasmaspheric  hiss  during  such  conditions  is 
plotted  as  a  function  of  L  shell  and  magnetic  local  time  in 
Figure  7.  The  hiss  intensities  peak  on  the  dayside  with 
intensities  typically  of  the  order  of  1000  pT2  over  a  range  of 
L  shells  from  L  =  2  to  L  =  4  from  0600  to  2 1 00  MLT.  Lower 
intensities  are  seen  near  midnight  and  at  larger  L  shells.  The 
latitudinal  variation  of  the  plasmaspheric  hiss  intensities 
averaged  over  all  MLT  are  shown  in  the  top  panels  of 
Figure  8.  The  hiss  intensities  tend  to  peak  near  the  equator 
(\m  <  5°)  and  at  higher  latitudes  (\m  >  15°). 

[32]  To  reduce  the  complexity  of  the  calculations,  we 
create  a  basic  model  of  the  plasmaspheric  hiss  intensities  as 
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Figure  7.  The  average  hiss  magnetic  field  wave  intensity 
as  a  function  of  L  and  MLT  for  quiet  conditions  ( Kp  <  3“). 
The  sampling  distribution,  color-coded  to  show  the  number 
of  minutes  in  each  bin,  is  shown  in  the  small  panel. 


a  function  of  L  shell  by  averaging  the  intensities  during 
quiet  conditions  in  steps  of  0.2 L  first  over  magnetic  latitude 
and  then  over  magnetic  local  time.  The  model  wave 
magnetic  field  intensities  are  plotted  as  a  function  of  L 
shell  in  Figure  9a  and  included  in  the  second  column  of 
Table  1.  The  model  intensities  have  a  peak  value  of 
1359  pT2  at  L  =  3.3  falling  to  1190  pT2  at  L  =  3.0  and 
380  pT2  at  L  =  5.0. 

[33]  The  mean  value  of  the  ratio  of fpJfce  is  plotted  as  a 
function  of  magnetic  latitude  for  the  same  conditions  and 
regions  in  the  bottom  panels  of  Figure  8.  Along  any  given  L 
shell  the  ratio  fpJfe  tends  to  decrease  with  increasing 


latitude  due  to  the  increasing  field  strength  with  magnetic 
latitude.  At  a  given  magnetic  latitude,  the  ratio  fpJfe  tends 
to  increase  with  increasing  L  shell.  The  dashed  lines  show 
the  model  values,  determined  from  least  squares  best  fits  to 
the  data,  assuming  a  dipole  field  and  constant  density.  The 
equatorial  values  of  the  best  fit  to  fpJfce  are  plotted  as  a 
function  of  L  shell  in  Figure  9b  and  are  tabulated  in  the  fifth 
column  of  Table  1.  The  equatorial  values  of  the  ratio  fpJfe 
increase  from  8.9  at  L  =  3.0  to  13.0  at  L  =  5.0. 

[34]  Plasmaspheric  hiss  appears  to  propagate  over  a  broad 
range  of  wave  normal  angles  with  predominantly  field- 
aligned  propagation  near  the  geomagnetic  equator  and 
more  oblique  propagation  at  higher  latitudes  [Parrot  and 
Lefeuvre,  1986;  Hayakawa  et  al.,  1986;  Santolik  et  al. , 
2001].  For  example,  in  the  equatorial  region  (\m  <  10°), 
Parrot  and  Lefeuvre  [1986]  found  two  populations  of  wave 
normal  angles,  one  lying  in  the  range  0°  <  \\)  <  30°,  the  other 
in  the  range  40°  <  \\)  <  60°.  At  higher  latitudes  (\m  >  20°), 
most  of  the  waves  had  larger  wave  normal  angles  in  the 
range  55°  <  o|j  <  85°.  To  investigate  the  effect  of  the  wave 
normal  angle  on  the  precipitation  lifetimes,  we  use  three 
different  angular  distributions  of  hiss,  chosen  to  be  repre¬ 
sentative  of  these  observations  (Table  2). 

[35]  The  conversion  from  electric  field  intensity  to  mag¬ 
netic  field  intensity  assumes  parallel  propagation  [Meredith 
et  al .,  2004].  We  calculate  approximate  intensities  for 
propagation  at  52°  and  80°  using  the  cold  plasma  dispersion 
solver  in  the  HOTRAY  code  [Horne,  1989]  assuming 
a  frequency  of  0.55  kHz  and  using  the  modeled  values  of 
fpJfce-  The  derived  intensities  are  tabulated  in  the  third  and 
fourth  columns  of  Table  1  and  are  typically  a  factor  of 
1.6  lower  for  a  wave  normal  angle  of  52°  and  a  factor  of 
5  lower  for  a  wave  normal  angle  of  80°. 

[36]  We  calculate  the  bounce-averaged  diffusion  rate 
which  takes  into  account  the  scattering  of  particles  in  pitch 
angle  over  the  complete  range  of  latitudes  between  the 
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Figure  8.  (top)  The  average  hiss  magnetic  field  wave  intensity  and  (bottom)  the  ratio  fpJf  e  for  quiet 
times  at  different  L  shells.  The  model  values  are  shown  by  the  dashed  lines. 
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Figure  9.  (a)  The  modeled  value  of  the  plasmaspheric  hiss  intensity  as  a  function  of  L  shell,  (b)  The 
modeled  value  of  the  ratio  fpJfce  at  the  magnetic  equator  as  a  function  of  L  shell. 


particle’s  mirror  points.  In  general  the  waves  resonate  with 
higher-energy  electrons  at  higher  latitudes  and  will  tend  to 
scatter  higher-energy  electrons  into  the  loss  cone  at  higher 
latitudes.  This  is  shown  in  more  detail  in  Figure  2  of  Horne 
and  Thorne  [2003]  for  the  case  of  chorus  waves.  We  assume 
a  Gaussian  frequency  distribution  centered  on  0.55  kHz 
with  a  bandwidth  of  0.3  kHz  and  lower  and  upper  cutoff’s 
of  0.1  kHz  and  2.0  kHz,  respectively. 

6.2.  Results 

[37]  Figure  10  shows  the  loss  timescale  as  a  function  of  L 
shell  for  214  keV  electrons  (lower  panel),  510  keV  electrons 
(middle  panel),  and  1 .09  MeV  electrons  (upper  panel).  The 
vertical  bars  denote  the  range  covered  by  plus  or  minus  one 
standard  deviation  of  the  mean  value  of  the  measured  loss 
timescales.  The  pluses  denote  the  model  values.  The  results 
are  shown  for  the  three  different  angular  distributions  of 
hiss,  blue  for  the  small  wave  normal  model,  green  for  the 
medium  wave  normal  model,  and  red  for  the  large  wave 
normal  model. 

[3s]  At  214  and  510  keV  the  small  and  medium  wave 
normal  models  produce  loss  timescales  that  are  consistent 
with  the  observed  values  in  the  regions  3.5  <  L  <  5.0  and 
3.0  <  L  <  5.0,  respectively.  At  1.09  MeV  both  models  fit 
the  data  well  in  the  region  3.0  <  L  <  4.0  but  overestimate 
the  loss  timescales  at  larger  L  by  a  factor  of  ~5.  In 
comparison,  the  large  wave  normal  model  overestimates 
the  loss  timescales  by  an  order  of  magnitude  or  more  at  all  L 
shells  and  energies  considered. 

[39]  The  lowest  loss  timescales  are  obtained  from  the 
small  wave  normal  model.  The  small  and  medium  wave 
normal  models  produce  loss  timescales  that  tend  to  lie 
within  a  factor  of  two  of  each  other,  whereas  the  large 
wave  normal  model  produces  loss  timescales  that  are 
typically  an  order  of  magnitude  larger. 

7.  Discussion 

[40]  The  energetic  outer  zone  electron  loss  timescales 
during  low  geomagnetic  activity  derived  from  the  CRRES 
MEA  data  compare  favorably  with  those  presented  in  previ¬ 
ous  studies.  For  example,  loss  timescales  derived  from  the 


OV3-3  satellite  in  the  region  3  <  L  <  4  following  a  period  of 
gradual  decay  following  the  intense  storm  of  4  September 
1966  range  from  2.6  to  5.1  days  for  475  keV  electrons  and 
from  5.1  to  8.5  days  for  957  keV  electrons  [Lyons  et  al. ,  1972]. 
Furthermore,  the  loss  timescales  derived  from  the  OGO-5 
satellite  in  the  region  3  <  L  <  5  following  the  moderate  storm 
of  11  June  1968  range  from  1.3  to  2.8  days  for  266  keV 
electrons,  from  1 .3  to  5.0  days  for  475  keVelectrons,  and  from 
1 .7  to  7.8  days  for  866  keV  electrons  [  West  et  al .,  1981].  More 
recently,  Baker  et  al.  [1994]  conducted  a  superposed  epoch 
analysis  of  1  year  of  SAMPEX  data  and  estimated  loss 
timescales  for  electrons  with  energies  greater  than  400  keV 
in  the  region  2.5  <  L  <  5  of  between  5  and  10  days.  While  a 
direct  comparison  with  a  specific  CRRES  energy  channel  is 
not  possible  in  this  case  the  range  of  values  has  some  overlap 
with  the  range  of  values  determined  from  CRRES  at  5 10  keV 
and  considerable  overlap  at  1.09  MeV  (Figure  10).  Albert 
[2000]  studied  the  9  October  1 990  storm  using  CRRES  MEA 
count  rate  data  and  showed  that  the  loss  timescale  of  5 10  keV 
electrons  had  a  minimum  of  1 .6  days  at  L  =  3.2,  rising  to  4.6 
days  at  L  =  2.6  and  2.6  days  at  L  =  3.7.  Seki  et  al.  [2005]  used 
CRRES  and  Akebono  data  to  estimate  the  loss  timescales  at 
L  =  3.25  following  the  strong  magnetic  storm  of  26  August 
1990.  The  loss  timescales  derived  from  the  CRRES  MEA  data 
ranged  from  to  3. 1  days  for  976  keV  electrons  to  5.3  days  for 
1.58  MeV  electrons  and  the  loss  timescale  derived  from  the 
Akebono  satellite  was  higher,  at  8.3  days  for  electrons  with 
energies  in  the  range  950-2500  keV.  The  longer  timescale 
observed  by  Akebono  was  attributed  to  the  slower  decay  of 
high-energy  electrons  within  the  instrumental  energy  range. 
The  lifetimes  derived  in  previous  studies  are  thus  largely 
consistent  with  the  range  of  lifetimes  determined  in  our  study 
(Figure  10),  both  as  a  function  of  energy  and  L  shell. 


Table  1.  Model  Parameters 
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Hadpl*) 
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3.0 

1190 
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10.0 
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4.0 
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Table  2.  Wave  Normal  Models 


Model 

xm 

hX 

*/, 

xm 

Small  wave  normal  model 

0.0 

0.00 

0.36 

0.00 

0.58 

Medium  wave  normal  model 

52.0 

1.28 

0.27 

0.84 

1.73 

Large  wave  normal  model 

80.0 

5.67 

2.74 

1.43 

11.4 

[41]  Albert  [1999]  calculated  plasmaspheric  lifetimes  for 
the  region  1 .5  <  L  <  5.0  using  the  parameter  values  of  Abel 
and  Thome  [1998]  for  whistler  mode  waves  due  to  hiss, 
lightning,  and  VLF  transmitters.  The  lifetimes  at  200  keV 
had  a  minimum  of  ~14  days  at  L  =  3.5,  increasing  to 
~50  days  at  L  =  3  and  ~20  days  at  L  =  5.  At  1  MeV  they 
had  a  minimum  of  ~50  days  at  L  =  3  rising  to  >400  days  at 
L- 2  and  ~200  days  at  L  =  5.  The  values  of  these  modeled 
lifetimes  are  typically  an  order  of  magnitude  larger  than  the 
experimental  values  derived  in  our  study.  For  L  >  3,  Abel 
and  Thorne  [1998]  used  a  hiss  amplitude  of  10  pT  assumed 
to  be  present  for  50%  of  the  time  corresponding  to  an 
average  intensity  of  25  pT2.  Our  measured  values  of  the  hiss 
intensities  have  a  peak  value  of  1359  pT2  at  L  =  3.3,  falling 
to  1190  pT2  at  L  =  3  and  380  pT2  at  L  =  5.  The  measured 
intensities  presented  here  are  thus  over  an  order  of  magni¬ 
tude  larger  than  the  hiss  intensity  used  by  Abel  and  Thorne 
[1998]  and  hence  also  by  Albert  [1999].  Since  lifetimes 
scale  inversely  with  wave  intensities,  the  reason  why  the 
plasmaspheric  lifetimes  obtained  by  Albert  [1999]  are  much 
larger  than  the  experimental  values  is  due  to  an  underesti¬ 
mation  of  the  hiss  wave  power. 

[42]  Our  results  suggest  that  plasmaspheric  hiss  propa¬ 
gating  at  small  or  intermediate  wave  normal  angles  is 
responsible  for  electron  loss  over  a  wide  range  of  L  shells 
and  energies  in  high-density  regions.  The  region  where  hiss 
dominates  loss  is  energy-dependent,  ranging  from  3.5  <  L 
<  5.0  at  214  keV  to  3.0  <  L  <  4.0  at  1.09  MeV.  The  results 
show  that  plasmaspheric  hiss  with  large  wave  normal  angles 
does  not  contribute  significantly  to  electron  loss  from  the 
radiation  belts  in  the  regions  shown. 

[43]  The  lowest  model  lifetime  for  214  keV  electrons  at 
L  =  3.3  is  a  factor  of  ~3  higher  than  the  experimental 
lifetime  at  L  =  3.35.  Unfortunately,  the  fitting  technique 
used  in  this  study  breaks  down  at  lower  L  for  214  keV 
electrons  due  to  increased  scatter  between  neighboring  data 
points  so  that  we  cannot  explore  the  discrepancy  further  in 
this  paper.  We  note  that  other  wave  emissions,  such  as 
lightning-generated  whistlers  and  VLF  transmitters,  become 
increasingly  important  at  low  L  and  will  contribute  to  the 
loss  timescales  in  this  region  [Abel  and  Thorne ,  1998]. 

[44]  Our  calculated  larger  loss  timescales  at  high  energies 
at  large  L  require  further  investigation.  Albert  [1994],  using 
a  highly  simplified  version  of  the  dispersion  relation  and 
resonance  condition  following  Lyons  et  al.  [1972],  found 
that  the  inclusion  of  higher-order  cyclotron  harmonics  could 
decrease  the  loss  timescales,  for  energies  above  500  keV 
and  L  values  greater  than  3.  In  particular,  he  found  that 
increasing  the  number  of  cyclotron  harmonics  from  ±5  to 
±100  decreased  the  1  MeV  loss  timescale  by  a  factor  of 
~2.3  at  L  =  5.  We  also  considered  the  contribution  from 
higher-order  cyclotron  harmonics,  increasing  the  number  of 
cyclotron  harmonics  to  ±100.  We  did  not  find  a  significant 
contribution  from  these  higher-order  cyclotron  harmonics 
for  any  of  the  models  used  in  this  study.  Higher-order 


cyclotron  resonant  effects  depend  sensitively  on  the  angular 
distribution  of  the  waves  and  the  ratio  fpJfce .  For  example, 
at  L  =  5  our  modeled  equatorial  value,  based  on  observa¬ 
tions,  for  the  ratio  fpJfce  is  12.8  compared  to  the  value  of 
26.9  used  by  Albert  [1994].  For  higher-order  cyclotron 
harmonics  (|/?|  >  5)  to  play  a  role  the  resonant  frequencies 
associated  with  these  harmonics  must  lie  within  the  fre¬ 
quency  band  of  the  waves.  In  our  models  the  wave 
frequencies  associated  with  these  higher-order  harmonics 
all  lie  above  the  upper  cutoff  frequency  for  the  range  of 
energies  and  L  shells  studied.  In  contrast,  in  the  study  of 
Albert  [1994]  the  wave  frequencies  associated  with  the 
higher-order  harmonics  fall  in  the  frequency  band  of  the 
waves  for  certain  energy  and  L  shell  ranges  and  hence  play 
a  role  in  the  diffusion  process.  Since  we  use  observed  values 
of  the  wave  parameters  and  the  ratio  fpJfce,  we  conclude  that 
higher  harmonic  effects  are  not  significant. 

[45]  Loss  by  chorus  waves,  known  as  microburst  precip¬ 
itation,  generally  occurs  outside  the  plasmapause  in  regions 
of  lower  density  [ Thorne  et  al .,  2005a,  2005b;  O  'Brien  et 
al .,  2004].  Thus  microburst  precipitation  could  contribute  to 
electron  loss  at  large  L  if  part  of  the  electron  drift  orbit 
occurs  outside  the  plasmapause.  However,  modeling  studies 
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Figure  10.  Electron  loss  timescales  as  a  function  of  L  shell 
for  (top)  E  =  1 .09  MeV,  (middle)  £=510  keV,  and  (bottom) 
£  =  214  keV.  The  measured  values  are  denoted  by  the 
vertical  lines  and  the  model  values  by  the  pluses,  color- 
coded  to  show  the  wave  normal  angle. 


11  of  13 


A05212 


MEREDITH  ET  AL.:  ENERGETIC  ELECTRON  LOSS  TIMESCALES 


A05212 


show  that  as  chorus  propagates  along  the  magnetic  field  the 
resonant  energies  match  ~  MeV  electrons  at  20° -30°  off 
the  equator  and  are  not  well  matched  at  the  equator  [Horne 
and  Thorne ,  2003].  At  the  equator  the  resonant  energies  are 
lower.  Thus  one  may  expect  that  if  microburst  precipitation 
were  occurring,  this  would  also  reduce  the  loss  timescales  at 
lower  energies.  It  appears  that  plasmaspheric  hiss  is  ade¬ 
quate  to  explain  the  loss  at  lower  energies.  This  does  not 
rule  out  microburst  loss  provided  the  rate  of  loss  at  lower 
energies  is  dominated  by  hiss. 

[46]  EMIC  waves  can  interact  with  MeV  electrons  [e.g.. 
Summers  and  Thorne ,  2003;  Albert,  2003;  Summers ,  2005], 
particularly  in  regions  of  high  density  at  larger  L  (L  >  4.5) 
[Meredith  et  al 2003c],  precisely  the  region  where  the 
MeV  losses  cannot  be  explained  by  resonant  interactions 
with  plasmaspheric  hiss  alone.  Furthermore,  Albert  [2003] 
found  that  the  addition  of  EMIC  waves  can  greatly  reduce 
the  electron  lifetimes  determined  by  hiss  alone.  In  addition, 
precipitation  signatures  from  EMIC  waves  have  been  seen 
in  plumes  during  substorm  ion  injections.  These  results 
suggest  that  EMIC  waves  in  high-density  regions,  such  as 
the  duskside  plasmapause  and  within  drainage  plumes,  may 
contribute  to  the  scattering  of  energetic  electrons  at  larger  L , 
and  could,  potentially,  explain  the  loss  rates  in  this  region. 

8.  Conclusions 

[47]  We  have  analyzed  CRRES  particle  data  to  determine 
the  energetic  electron  decay  timescales  as  a  function  of 
energy  (2 14  keV  <  E  <  1 .09  MeV)  and  L  shell  (3  <  L  <  5). 
The  resulting  timescales  are  compared  with  estimates  of  the 
decay  timescales  of  plasmaspheric  hiss  using  the  PADIE 
code  using  CRRES  wave  data  as  input.  Our  principal  results 
are  as  follows: 

[48]  1 .  Gradual  loss  of  energetic  electrons  in  the  region 
3  <  L  <  5  is  observed  to  occur  in  quiet  periods  (Kp  <  3~) 
following  enhanced  magnetic  activity. 

[49]  2.  Decay  timescales  range  from  1.5-3. 5  days  for 
214  keV  electrons  to  5. 5-6. 5  days  for  1.09  MeV  electrons. 

[50]  3.  The  intervals  of  decay  are  associated  with  large 
average  values  of  the  ratio  fpJfce  (>7)  and  indicate  that  the 
decay  takes  place  in  the  plasmasphere. 

[51]  4.  Modeling  shows  that  plasmaspheric  hiss  propa¬ 
gating  at  small  or  intermediate  wave  normal  angles  can 
account  for  the  electron  loss  rates  over  a  wide  range  of 
energies  and  L  shells.  These  regions  are  energy-dependent 
and  range  from  3.5  <  L  <  5.0  at  214  keV  to  3.0  <  L  <  4.0 
at  1.09  MeV 

[52]  5.  Modeling  shows  that  plasmaspheric  hiss  with  very 
large  wave  normal  angles  does  not  contribute  significantly 
to  radiation  belt  electron  loss. 

[53]  6.  At  high  energies,  E  =  1 .09  MeV,  plasmaspheric  hiss 
propagating  at  small  or  intermediate  wave  normal  angles 
overestimates  the  loss  timescales  by  a  factor  of  ~5  at  large  L 
(4.5  <  L  <  5.0).  Resonant  wave-particle  interactions  with 
EMIC  waves  become  important  at  MeV  energies  at  larger  L 
and  high  densities.  We  suggest  that  losses  due  to  EMIC  waves 
may  play  a  significant  role  at  larger  L. 
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